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Impact of Dynamical Fermions on QCD Vacuum Structure 
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We examine how dynamical fermions affect both the UV and infrared structure of the QCD 
vacuum. We consider large 28^ x 96 lattices from the MILC collaboration, using a gluonic definition 
of the topological charge density, founded on a new over-improved stout-link smearing algorithm. 
The algorithm reproduces established results from the overlap formalism and is designed to preserve 
nontrivial topological objects including instantons. At short distances we focus on the topological 
charge correlator, {q{x)q{0)), where negative values at small x reveal a sign-alternating layered 
structure to the topological-charge density of the QCD vacuum. We find that the magnitudes of the 
negative dip in the {q{x)q{0)) correlator and the positive {5(0)^) contact term are both increased 
with the introduction of dynamical fermion degrees of freedom. This is in accord with expectations 
based on charge renormalization and the vanishing of the topological susceptibility in the chiral 
limit. At large distances we examine the extent to which instanton-like objects are found on the 
lattice, and how their distributions vary between quenched and dynamical gauge fields. We show 
that dynamical gauge fields contain more instanton-like objects with an average size greater than 
in the quenched vacuum. Finally, we directly visualize the topological charge density in order to 
investigate the effects of dynamical sea-quark degrees of freedom on topology. 

PACS numbers: 12.38.Gc ll.15.Ha 12.38.Aw 



I. INTRODUCTION 

The study of QCD vacuum structure is one area of 
research where lattice simulations provide access to oth- 
erwise unaccessible information. By generating typical 
vacuum gauge field configurations we are able to directly 
investigate their complex structure. Further, by varying 
the simulation parameters we can assess how different 
physical phenomena contribute to the structure of the 
vacuum. In particular, we are able to examine the effect 
of dynamical sea-quarks on the QCD vacuum. 

A study of QCD vacuum structure at different scales 
usually requires the use of a filtering procedure. By re- 
moving the short-range UV fiuctuations, one can probe 
the long-distance structural features of the vacuum. 
Typical UV filtering methods include coohng P, 0, [1], 
APEQI, I^] , and improved APE smearing @ , HYP smear- 
ing [31 and more recently, stout-link smearing Q and 
LOG smearing |9|. More recently, the truncation of a 
spectral representation of operators using eigenmodes 
of the Overlap Dirac operator has been explored as a 
method to remove UV fluctuations [H [H [ll [11] . 

Instantons are believed to be an essential component of 
the long-distance physics of the QCD vacuum. The prob- 
lem with most cooling and smearing algorithms is that 
they destroy the instantons in the vacuum. In particular, 
standard EXP, or stout-link, smearing suffers from this 
problem. In order to overcome this we employ a form 
of over- improved stout-link smearing [l5| , described 
briefly in Sec. |TT1 

One facet of recent topological studies of vacuum struc- 
ture is the topological charge density correlator, or Eu- 
clidean 2-point correlation function 



which is the integrand of the topological susceptibility 



V 



= d^x{qix)q{0)), 



(2) 



where V is the 4-volume. 

In Sec. mil we examine the shapes of the (qq) corre- 
lator, using both over-improved stout-link smearing and 
3- loop, 0{a'^) improved cooling. We use the gluonic def- 
inition of the topological charge density 



q{x) 



9 



where 



(3) 



(4) 



{qq) ^ {q{x)qm , 



(1) 



and Ff^ij (x) is a 3-loop improved field strength tensor . 
Of particular interest is the extent to which a purely glu- 
onic definition of the topological charge density can revel 
a negative (qq) correlator. We then proceed to inves- 
tigate how the correlator is affected by the inclusion of 
dynamical sea-quarks. 

In Sec. IIVI we probe the infrared structure of the full 
dynamical QCD vacuum. We examine the extent to 
which instantons are found in the QCD vacuum and how 
the extra dynamical fermion degrees of freedom affect 
their distribution throughout the vacuum. 

Finally, in Sec.|V]we directly visualize the topological 
charge density in order to investigate the effects of dy- 
namical sea-quark degrees of freedom on topology. Con- 
clusions are summarized in Sec. IVIl 
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II. OVER-IMPROVED STOUT-LINK 
SMEARING 

Extended smoothing on any QCD gauge field will re- 
veal the presence of smooth topological objects that are 
approximations to the classical Euclidean instanton so- 
lution. A single instanton has the gauge potential [iTj 



A,{x)^ , (-]d,{S)S-\ (5) 



where 



S = 



Xi ± i X ■ a 



(6) 



In order to tame the lattice discretisation errors, 
improved actions have been developed, such as the 
Symanzik [2ll | 0(a^ )-improved action 



12 



((l-i?^,(a;)) + (l-i?,^(a;))) 



(10) 



where we have and to denote the different possi- 
ble orientations of the rectangular loop, and P^^ denotes 
the plaquette. 

We modify this action, following Perez et al. [lO], by 
the inclusion of a new e parameter 



for instantons and anti-instantons. 

Due to their correlations with low-lying Dirac eigen- 
modes these locally self-dual instanton-like objects are of 
obvious physical interest. Unfortunately it is well real- 
ized that standard smoothing procedures will often de- 
stroy instantons in the gauge field. 

For this investigation we have chosen to employ over- 
improved stout-link smearing ^14]. and will now present 
a brief summary of the algorithm. Complete details 
and an investigation of the algorithm's performance are 
provided in Ref. [l3]- We also present results using 3- 
loop iinproved cooling [l^, and defer the reader else- 
where for a discussion of highly improved 
actions. 

Standard stout-link smearing is closely related to the 
Wilson action 



Sw{x) = 13 



E 

1/ 



^RetT[l-U^ix)Ef,,ix)] , (7) 



where 



^f^Ax) = U,ix + A) Ulix + v) Ulix) 

+ Ul{x + (i-v) UUx - v) U^x - v) , 



(8) 



to calculate the nearest neighbor contributions. However, 
it is possible to use any reasonable combination of links 
connected to the original link [/^t (x) . We can analyze the 
effect of the Wilson action on an instanton by substitut- 
ing the instanton solution into the Taylor expansion of 

Sw{x) [a [3 



9 



1 



5 VPinst/ 70 VPinst 



(9) 



If one plots this expression as a function of the instanton 
size pinst one sees that the slope of the action is necessar- 
ily always positive. Denoting the action associated with 
a single classical instanton as Sq, we note that the ideal 
action for preserving instantons in the gauge fields would 
have a flat curve at 5/50 = 1, as in the continuum. 



^w = /3EE 

X l-L>V 

1 - e 



2e 



12 



(11) 



where for e = 1 we have the standard Wilson action, and 
for e = we have the Symanzik improved action. So by 
varying e between and 1 we can effectively select how 
much improvement is included. 

Alternatively we can set e < and over-improve the 
action. The advantage of a negative e value is that the 
leading order errors become positive upon substitution of 
the instanton solution 



9^ 



e / a 
1 - - ' 



5 VPinst 



14e- 17 / a 



210 VPinst 



(12) 

The curve of S/Sq is similarly affected. One sees that 
for a value of e = —0.25 the curve is mostly flat and 
will preserve instantons objects on the lattice with size 
Pinst/a > 1.5 [H]. 

Standard stout-link smearing, using an isotropic 
smearing parameter psm, involves a simultaneous update 
of all links on the lattice. Each link is replaced by a 
smeared link ?7p(x) 



where 



with 



flf_,{x) 



Uf,ix) = exp(iQ^(x)) U^{x) , (13) 



QAx)^-{nl{x)-n^{x)) 

-'-Tr{nl{x)-n,{x)), (14) 



y~^{l X 1 loops involving U^{x)} . (15) 



The over-improvement parameter e is introduced into the 
smearing process by replacing the combination of links in 



3 



flfj,{x) with 



^ f 5 — 2e 

= psm ^ I — ^ — (1 X 1 loops involving 



1 - e 



1x2 + 2x1 loops involving U^{x)) 



(16) 



Note that both forward and backward horizontally orien- 
tated rectangles are included in the 2x1 loops, such that 
ri^(x) resembles the local action. With this extended 
link path, we take — 0.06, smaller than the standard 
value of 0.1. 



III. TOPOLOGICAL CHARGE DENSITY 
CORRELATOR 

A. Quenched QCD 

The Euclidean 2-point correlation function for the 
topological charge q{x) 



{qq) ^ {q{x)qm , 



(17) 



also known as the topological charge density correlator, 
is negative for any x > [23 , . This follows simply 
from reflection positivity [24| . Given that the correlator 
must have a positive contact term {q'^{0))x, the correlator 
necessarily has the form 



{q{x)qm = AS{x) - fix) , 



(18) 



where f{x) is positive for all x. 

Recent studies [lol . [24| of vacuum structure using the 
overlap defined topological charge density have demon- 
strated the negativity of the topological charge density 
correlator in lattice simulations. The negative behavior 
of the (qq) correlator suggests a sign-alternating layered 
structure to the topological charge density of the topo- 
logical charge density correlator. 

Hasenfratz [2^ has also examined the shape of the 2- 
point correlator using APE smearing. For APE smear- 
ing the magnitude of the negative dip was far less than 
the current overlap results [lOl . [23 |. We now attempt 
to obtain a negative (qq) correlator similar to the over- 
lap results using the new over-improved stout-link smear- 
ing algorithm and 3-loop improved gluonic operators for 
the action and topological charge densities . In order 
to more closely compare with the overlap results we use 
quenched gauge fields for the initial calculation. 

The gauge fields were generated by the MILC collabo- 
ration [2a, 123] using a Tadpole and Symanzik improved 
gauge action with 1x1-1-1x2-1-1x1x1 terms in 
the quenched case and an Asqtad staggered dynamical 
fermion action for the 2 + 1 fiavors of dynamical quarks. 
The lattice spacing for all three types of gauge fields is 
a — 0.09 fm. For the specifics of how the gauge fields 



TABLE I: The gauge fields used in this st udy . The lattices 
were generated by the MILC collaboration [2a . [27j 



Size 


/3 


a 


Bare Quark Masses 


28^ X 96 


8.40 


0.086 fm 


oo 


28^ X 96 


7.09 


0.086 fm 


14.0 MeV, 67.8 MeV 
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0.086 fm 


27.1 MeV, 67.8 MeV 
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FIG. 1: The {qq) correlator for a variable number of smooth- 
ing sweeps calculated on the quenched gauge fields. The top 
graph shows the (qq) correlator for over-improved stout-link 
smearing smearing, and the lower graph for 3-loop improved 
cooling. Using four sweeps of smearing and one sweep of cool- 
ing it is possible to obtain a correlator similar to the recent 
overlap results [l^, [l^] . As more U V fluctuations are removed 
in successive smoothing iterations the correlator flattens off. 
Errors were calculated using the jackknife procedure but are 
too small to see. 



were generated see Refs. [26|,[27[. Simulation parameters 
are summarized in Table [Ij 

We begin by investigating how the shape of the corre- 
lator changes as a function of the number of smearing 
sweeps. The results of our calculations are presented 
in Fig. [1] We see that for a small number of sweeps 
it is possible to generate a negative correlator, but that 
the negative behavior is largely suppressed after about 
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FIG. 2: The topological charge density of a quenched gauge 
field after four sweeps of smearing. Note the sign-alternating 
structure in the charge density. 



10 sweeps and absent after around 40 sweeps. 

The curves also have a positive core whose minimal size 
is limited by the lattice spacing and the combination of 
loops in the smoothing procedure. As more of the short- 
distance UV fluctuations are removed by the smoothing 
algorithm there is a dampening effect on the topological 
correlator. 

One sweep of highly improved cooling gives a (qq) cor- 
relator of similar proportions to the overlap defined cor- 
relator , whilst for over-improved smearing, four 
sweeps reveals a similar curve. A quantitative compar- 
ison of the overlap correlator is performed in Ref. 
The topological charge density of a quenched field after 
4 sweeps of smearing is presented in Fig. [2] There is 
considerable structure within the gauge field. A sign- 
alternating structure in the topological charge density is 
also apparent. 



Including dynamical sea-quarks in the QCD action 
renormalizes the coupling constant. In order to maintain 
the same lattice spacing across quenched and dynami- 
cal gauge fields, one finds that the coupHng parameter, 
g, must increase and hence (3 ~ 1/5^ must be smaller 
for the dynamical fields. In the QCD action formulated 
in Euclidean space, (3 appears as a factor governing the 
width of the probability distribution. When generating 
gauge fields, smaller (3 values permit greater fluctuations 
in the gauge links. The increased fluctuations can give 
rise to non-trivial field fluctuations, which will reveal it- 
self in lattice simulations through a greater mean-square 
density {q'^iO))^. 

Given that in the chiral limit )(^y" = 0, it follows 
that an increasing mean-square density {q'^{0))x must 
be compensated for by a stronger negative dip in the 
(qq) correlator. The (qq) correlator for the quenched 
ensemble and the two dynamical ensembles is shown 
in Fig. [31 These correlators were generated using four 
sweeps of over-improved stout-link smearing. We see that 
the contact term (^^(O))^; is larger and the magnitude of 
the negative dip has also increased. These effects are 
stronger for lighter quark masses. This is of course ex- 
pected because smaller quark masses require smaller /3 
values. The exact values of the positive contact term 
are; quenched = 2924 ± 4 fm"^, heavy = 5251±12fm"*, 
light = 5432 ± 8 fm'^ 

Hasenfratz [2^ also observed an increase in the nega- 
tive contribution to x with the addition of fermion loops. 
However, their correlators also show a decrease in the 
mean square density, which is in contrast to our results. 
This may be explained through the use of a standard 
smoothing algorithm, which is renowned for destroying 
topological objects in the gauge field. 



IV. INSTANTON-LIKE OBJECTS 



B. QCD with Dynamical Fermions 



A. Profile Versus Charge Density 



We now investigate the effects of dynamical sea-quark 
degrees of freedom on topology. In quenched QCD the 
Witten-Veneziano [H, [2§| formula gives a relation be- 
tween the topological susceptibility and the mass of the 
ij' meson j23| 



X 



2Nf 



(19) 



However, in the full dynamical case the topological sus- 
ceptibility should vanish in the chiral limit [23l. |30|| 



dyn 



2Nf 



(20) 



Of course, a vanishing topological susceptibility puts 
no restraints on how the function {q{x)q{Q)) should 
change with the addition of dynamical sea-quarks, it only 
requires that the integral of Eq. ^ vanish. 



Understanding the nature of instanton-like objects in 
the QCD vacuum continues to be an active of area of 
investigation. Considerable UV filtering reveals the pres- 
ence of long-distance topological structures in the QCD 
vacuum. These topological objects are only approxima- 
tions to the classical instanton solution, but are com- 
monly referred to as instantons. The amount of smooth- 
ing usually applied in order to observe these instanton- 
like objects is far greater than the four sweeps of over- 
improved stout-link smearing that we used in the calcu- 
lation of the {qq) correlator. 

We now analyze the similarity of the topological ob- 
jects in the QCD vacuum to the classical instanton so- 
lution. Using over-improved stout-link smearing we are 
able to extract both the action and charge densities of 
our gauge fields. Starting with the action density we lo- 
cate the positions of all local maxima in the field. The 
local maxima are identified by finding a point at the cen- 
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FIG. 3: A comparison of the (qq) correlator for quenched and 
dynamical gauge fields. The smaller /3 values used for the 
dynamical fields result in greater field fluctuations that are 
visible through an increase in the magnitude of the negative 
dip. Although not shown, the contact term {q^{0))x has also 
increased and values are given in the text. The lower graph 
displays the same data as the upper graph, focusing on the 
negative dip for the dynamical fields. The magnitude of the 
dip is greater for the lighter quark mass. 



ter of a 3'* hypercube whose action density exceeds that 
of the neighboring 80 points of the hypercube. 

Taking each maxima to be the approximate center of a 
possible instanton-likc object we fit the classical instan- 
ton action density 



6 



Pfnst 



{{x - xoY 



Phist ) 



(21) 



to the measured action density. An arbitrary scale factor, 
^, is included to allow the shape of the action density to 
determine the size, Pinst- ^Ve fit the six parameters, xg, ^ 
and pinst by fitting Eq. (|21|) to the action density of the 
aforementioned 3* hypercube. 

From pinstj one can infer the topological charge to be 
observed at the center of the distribution q{xo) if it truly 



q{xo) = Q , 4 

^ Pfnst 



(22) 



Here Q = ±1 for an instanton/anti-instanton. This can 
then be compared with the topological charge measured 
directly from the charge density observed on the lattice. 

A calculation on the ensemble of quenched gauge fields 
is presented in Fig. H) Each cross represents a peak in the 
action density. If the peak were to exactly represent an 
instanton it would lie on the theoretical curve of Eq. (|22|1 . 
For a small number of sweeps we see that the peaks in 
the action density do not correspond well to instantons. 
It is only after extensive smearing has been applied that 
the objects start to approximate instantons. 

After 40 sweeps of over-improved stout-link smearing 
we have a large number of maxima in the vacuum which 
appear to be good approximations to the classical in- 
stanton solution, at least within the S** hypercube used 
to determine the "instanton" size. 

While the core shape and density approximate an in- 
stanton, we now assess how similar the remainder of 
the object is to a classical instanton. To do so, we 
consider all points within a distance r measured rela- 
tive to the instanton size pinst and examine the extent 
to which the points within this distance have the same 
sign topological-charge density as observed at the center, 
q{xo). If the detected object is a good approximation to 
a classical instanton then all these points should have the 
same charge as q{xo). 

In Fig. [5] we show the percentage of points that are 
within a relative distance r/pinst of xq that have the 
same sign as xq. For small r/pinst the percentage of sign- 
coherent points is close to 100, however it falls off rapidly 
as r approaches the characteristic size pinst • This suggests 
that although the object is representative of an instan- 
ton at its center, the tails of the objects are distorted 
by vacuum fluctuations. What is remarkable is that at 
the characteristic size of the "instanton," merely 2/3 of 
the points are sign coherent, suggesting that the objects 
revealed after 40 sweeps of smearing are good approxi- 
mations of classical instantons only at the core. 



B. Dynamical Fermions and Instanton Structure 

Given the strong correlation between q{xo) extracted 
from the action density after 40 sweeps of smearing and 
that given by the classical solution we can now com- 
pare "instanton" distributions between quenched and dy- 
namical QCD. 

We examine the variation in instanton size between 
the different gauge fields. A histogram of pinst is shown 
in Fig. [51 Compared to the quenched ensemble, the dy- 
namical gauge fields show an increase in average instan- 
ton size and also in the total number of instantons. A 
similar comparison was carried out in Ref. [3l| where 
it was argued that the instanton distributions were the 
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FIG. 4; A comparison of the charge at the center of the 
instanton-hke objects q(xQ) observed on the lattice versus 
their size pinat obtained via Eq. (I2ip for three different levels 
of over- improved stout-link smearing. From top to bottom 
we apply 4 sweeps, 20 sweeps and 40 sweeps of smearing. 
Each cross represents a local maxima in the action density of 
the gauge field. After 40 sweeps most points lie close to the 
classical curve. 



same between quenched and dynamical gauge fields, how- 
ever the results of that study were limited by statistical 
fiuctuations. 

It is often argued that the increased density of 
instanton-like objects on the lattice can be explained 
through an instanton/anti-instanton attraction occurring 
due to the presence of the fermion determinant in the 
QCD weight factor [sO]. Isolated instantons and anti- 
instantons give rise to zero modes of the Dirac operator. 



^ 60 




FIG. 5: The percentage of points that are sign coherent within 
a relative distance r j p^-ast of each instanton-like objects' center 
xq with size pinst. For small r j pi^^x. the percentage of sign- 
coherent points is close to 100, however the sign-coherence 
falls off rapidly as r approaches the characteristic size pinat. 
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FIG. 6: A histogram of the instanton size pinat in dynamical 
and quenched gauge fields. The dynamical fields show both an 
increase in the number of instantons detected and in average 
instanton size. 



When generating dynamical gauge fields the selection of 
typical configurations is weig hted by det(0 + 
If an exact zero-mode of ij) were to exist on the lattice 
then the determinant would approach in the chiral limit 
and it would be highly improbable that the configuration 
would be selected. Thus, isolated instantons will not ex- 
ist in the light quark mass gauge fields and hence all 
instanton-like objects will be closer in these fields. This 
naturally leads one to expect that there will be a greater 
number of instantons in the dynamical gauge fields. How- 
ever, this docs not account for the increased average size 
of the topological objects, which is an interesting result. 
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V. TOPOLOGICAL CHARGE DENSITY 

The effects of dynamical quarks are realized signifi- 
cantly in the calculation of the (qq) correlator. The in- 
creased magnitude of the non-trivial topological charge 
field fluctuations that are permitted due to the inclusion 
of fermion loops should also be visible in direct visualiza- 
tions of the topological charge density. 

Using four sweeps of over-improved stout-link smear- 
ing, we consider the short-range structure of the topolog- 
ical charge density. In Fig. [7| we present the topological 
charge density for the quenched and two dynamical en- 
sembles. The extra field fluctuations that are permitted 
due to the fermion determinant are clearly visible in the 
visualizations of the dynamical QCD vacuum. 

In Fig. [5] we compare the structure of the vacuum after 
40 sweeps of over-improved smearing. It is difficult to ob- 
serve the increased density of instantons in these figures. 
This is because the charge density fiuctuates over dif- 
ferent time-slices and a single time-slice does not always 
portray an accurate representation of the vacuum. 

VI. CONCLUSION 

By applying two different levels of smearing to the 
large-volume MILC gauge fields, we explored the effects 
of dynamical fermions on topology. The results agree 
with expectations. The addition of fermions into the 
QCD action at constant lattice spacing renormalizes the 
coupling constant such that the coupling parameter /3 
becomes smaller. This permits greater field fiuctuations 
and is manifest through an increased RMS density of 



topologically non-trivial field fluctuations. This induces 
an increased negative dip in the topological charge den- 
sity correlator and a larger contact term reflecting the 
RMS density of topological charge density. This effect 
increases as rriu^d ^ 0. 

The results reflect the suppression of zero-modes due 
to the inclusion of the det(0 -I- m) weight factor in the 
selection of typical gauge fields, resulting in a decrease in 
the number of isolated instanton-like objects. This causes 
instantons and anti-instantons to be "attracted" [3(?| and 
leads to an increase in the density of instanton-like ob- 
jects in the dynamical gauge fields. 

These results support the emerging picture of the vac- 
uum as an alternating "sandwich" of opposite topo- 
logical charge density. Beneath this oscillating short- 
range structure there exists a long-distance foundation 
of instanton-like objects that can be revealed through 
smoothing. The addition of dynamical fermions allows 
stronger field fluctuations and a higher frequency of sign- 
oscillations in the topological charge density. The density 
of instanton-like objects beneath these short-distance os- 
cillations also increases, as does their average size. 
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FIG. 8: The topological charge density for the quenched and dynamical ensembles, obtained after forty sweeps of over-improved 
stout-link smearing. From top to bottom we plot a quenched field, the dynamical am^^d = 0.0124, arris = 0.031 field and the 
dynamical amu,d = 0.0062, arris = 0.031 field. It is difficult to observe a noticeable difference in the size-distribution of the 
gauge field fluctuations when observing a single time-slice, but the sparseness of the quenched field is readily observed. 



